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Abstract — A general theory of inhomogeneous broadening is 
rarely applied to nonlinear spectroscopy in lieu of either a simple 
Lorentzian or Gaussian model. In this work, we generalize all 
the important third-order nonlinear susceptibility expressions 
obtained with sum-over state quantum calculations to include 
Gaussian and stretched Gaussian distributions of Lorentzians. 
This theory gives a better fit to subtle spectral features - 
such as the shoulder of the electroabsorption peak, and is a 
more accurate tool for determining transition moments from 
spectroscopy experiments. 



I. Introduction 

SUM-over states quantum perturbation treatments of the 
b th -order nonlinear susceptibility tensor, £y k , in the 
dipole approximation yields a sum of terms of the form: 
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where (n>i) nm is me nm-matrix element of the i th component 
of the electric dipole operator, uj nm the transition frequency 
(energy) between states n and m, Ui the frequency of the i th 
optical field, and r„ 9 the phenomenological damping factor. 
The numerator is a product of b + 1 transition moments and 
the denominator a product of b energy terms. For an isolated 
molecule, the damping factor T ng is inversely proportional to 
the lifetime of state n and is a measure of the width of the peak 
in the spectrum of k associated with a transition between 
state n and the ground state g. 

In this work, we take into account the interaction 
of molecules with their surroundings using a stochastic 
model as we have reported for the first- and second-order 
susceptibility.] 1 1 The technique is similar to Stoneham's ap- 
proach for the linear susceptibility |2| and Toussaere's calcu- 
lation of the hyperpolarizability,[3| who both used Guassian 
statistics. In our work, we generalize the statistics to stretched 
exponentials, which are known to better model the interaction 
between a molecule and a system that is characterized by a 
distribution of sites such as a host polymer. [4], [5| 

In our treatment of inhomogeneous broadening, each 
molecule in an ensemble is then viewed as having a different 
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transition frequency (energy), u> ng - For the stretched Gaussian, 
the probability distribution is of the form 
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transition frequency, N(j ng ,(3) the normalization factor, 7„ 9 
the linewidth of the distribution and (3 is the distribution of 
sites parameter. For most systems, (3 = for an infinitely 
broad distribution and (3 = 1 for a single characteristic width. 
The susceptibility will then be of the form, 
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Note that AT (7, (3) depends on j3, and will be written as 

N( lng ,(3)=~f ng V^B((3), (4) 
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to remain compatible with previous inhomogeneous broaden- 
ing representations that use Gaussian statistics [3|, [6|, [7|, [8|. 
Note that such stretched Gaussian statistics are also observed 
in light scattering experiments[8j, (9). 

In this paper, we derive the expressions for the most 
important third-order susceptibilities. The specific results for 
a Gaussian and stretched exponential are presented and com- 
pared with data from a quadratic electroabsorption (third-order 
susceptibility) experiment to illustrate the usefulness of the 
theory. All results are summarized in an extensive appendix. 

II. Third-Order Energy Denominators 

Similar to first- and second-order processes! 1], the SOS 
Lorentzian energy denominators for third-order processes are 
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In our notation, Df n {— u>i, 0^2, ^3) represents 
interactions, which involve only one-photon states, and 
(— Wi, W2, k/3) represents interactions that involve 
both one- and two-photon states. 

It is significantly more difficult to calculate the third-order 
IB theory because of the triple product of Lorentzian terms 
in the denominator. In order to transform Df , and/or Df 
for a specific set of input and output frequencies, the number 
of excited states must be known prior to performing a partial 
fraction expansion of each denominator term. For example if 
there are two distinct excited states (I and n), it is necessary 
to perform the following partial fraction expansion, 
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to eliminate the product of the two fli g terms. These type of 
expansions allow us to write the nonlinear energy denomina- 
tors in terms of wi 1 ^ (z) or when /3 = 1 in terms of complex 
error functions. 

The first and second-order IB theory has been derived from 
the Lorentzian denominator terms by Rruhlak and Kuzyk,[l| 
so we begin with the analogous derivation of the fundamental 
transformation for the third-order denominator to model all 
third-order processes. As an example consider, 
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Equation (O is multiplied by the stretched Gaussian function 
and is integrated with respect to 8u> ng , the integration variable 
is changed to t = (cj' ng - uj ng )/j ng , and z replaces (—cj ng + 
i^ng + ^)/lng, to get the following: 
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where 
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and Z = -{(jJng ~ iT ng ~ u)Hng< 

The integral in Eq. ( fTTT ) looks similar to the complex error 
function] 10), 

exp(-i 2 ) 
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when (3—1, except for the denominator. We use integration 
by parts to re-express the denominator to first-order in (z — t). 
With T = exp(— t 2 ), integrating by parts twice yields: 



T 



rdt 



_^M. {Z-tf J_^JL (Z-t) 



tT 



rdt 



= 2 



t 2 T 



-dt 



"ng ( Z - t) J^^UL (Z - t) 



T 



-dt, (13) 



where certain terms vanish when the argument of the expo- 
nential is small (w — 10 3 ) at the lower limit. 

Using (z+t) = (z 2 -t 2 )/(z-t) and t exp(-t 2 )dt = 
to recast Equation ( fT3l ) into a more convenient form, we get: 
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So, the convolution of the cubic Lorentzian with the Gaussian 
distribution (with (3 = 1) is, 
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where z = (-0J ng + i^ng + v)/jng- 

Table J] summarizes the third-order fundamental energy 
denominators for the Lorentzian and IB theories with (3 < 1, 
and (3—1, respectively. These in conjunction with those 
derived in Ref. [1J can be used to construct any IB energy 
denominator for any first-, second-, and/or third-order process. 



III. Third-Order Molecular Susceptibility 

As an example of the third-order molecular susceptibility 
for homogeneous and inhomogeneous models, we use Eq.'s 
©, Q, (ISgll,®,® and ([32}, for the respective models, in 
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to calculate the imaginary part of £( 3 ) (— u>; to, 0, 0) for a three- 
level system. Figure Q] shows the imaginary part of the third- 
order susceptibility for a system with one one-photon and 
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TABLE I 

Fundamental denominator contributions to homogeneously 
broadened and inhomogeneously broadened electronic 
transitions 
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Fig. 1. Imaginary part of £( 3 ) (— u; £j, 0, 0) from the generalized IB theory 
for a one-photon excited state centered at 660 nm (Ti g = 10 meV, and 
■yig = 40 meV) and a two-photon state centered at 595 nm (T2 g = 40 meV, 
and 729 = 40 meV). Three values of are compared to the Lorentzian theory 
(ri 9 = 40 meV, and F2 g = 40 meV). H2 g /fJ,i g = 0.4 for all the models. 



one two-photon excited state for the quadratic electrooptic 
process. We see a decrease in the magnitude of the third-order 
susceptibility from the Lorentzian theory to the IB theory, 
which suggests that inhomogeneous broadening decreases the 
effective nonlinearity of the material for this process. This 
trend continues as the distribution of sites becomes broader. 
Additionally the significant reduction in the response near the 
center wavelength of the two-photon state can be very dramatic 
when the contributions from the one- and two-photon states 
are of similar strength. This occurs because they have opposite 
signs in Eq. ( fl6b . Therefore it is important to use the IB theory 
for the guest-host system to model this particular region. 

IV. Comparison of Theory to Quadratic 
Electroabsorption Experimental Results 

We use quadratic electroabsorption spectra to test our mod- 
els. Details of the experiment and the relationship between 
Eq. (O and x (3) ( — us;us, 0,0) can be found in the literature. 
ifTTl Figure [2] compares experimental values of the imaginary 
part of for silicon phthalocyanine-methylmethacrylate in 
polymethylmethacrylate (SiPc/PMMA), which were derived 
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Fig. 2. Absolute value of the imaginary part of x' 3 ' from a quadratic 
electroabsorption experiment on SiPc/PMMA and least-squares fits using 
Lorentzian and IB(/3 = 1) theories. The fit parameters are from ref. 1 1 1 1. 



from quadratic electroabsorption experiments. Also plotted are 
the Lorentzian (L) and IB theories. A log scale is used to 
highlight the qualitative and quantitative features of the two 
models. The fit parameters are from the literature] 11 1. Like 
the linear absorption fits in Ref. (TJ, the IB model fits the data 
better both quantitatively, roughly a factor of 6 smaller relative 
error, and qualitatively, especially in the wings of the resonant 
signal. Neither model fits the data off-resonance because of 
the large random error associated with the lock-in amplifier 
signal away from resonance. The error bars cannot be plotted 
on a log scale because the error range includes negative values 
in the wings. 

Nonlinear spectroscopy experiments aim to determine zero 
frequency nonlinear susceptibilities by extrapolation, which 
can lead to large uncertainties depending on the quality of 
the dispersion models. Indeed, Canfield lTl2l . Vigil. fOl and 
Kruhlak IfTTl have shown that it is often difficult to reconcile 
the transition moments as determined by independent means. 
The wing region near resonance and the shape of the resonance 
peak may play an important role when using fitting to deter- 
mine transition moments or for extrapolating to off-resonant 
values of x from a data set with limited spectral range. 

Differences between the IB model and the standard 
Lorentzian model can be used to determine the reliability of 
zero-frequency susceptibilities and the uncertainty in transition 
moments. More importantly, precise modeling aimed at under- 
standing the dispersion of the nonlinear-optical response will 
need to take into account all possible broadening mechanisms. 
Because IB theory takes into account the distribution of sites, 
it may well be the best way to model systems such as dye- 
doped polymers. 

V. Conclusion 

In conclusion, we have calculated the inhomogeneously 
broadened third-order nonlinear-optical susceptibilities for a 
Gaussian and stretched Gaussian distribution of Lorentzians. 
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The results are applied to the quadratic electro-absorption 
spectrum of SiPc/PMMA and we find that the Lorentzian 
fit alone does not fit the data at the wings. The IB theory, 
however, fits the data over a broader wavelength range and 
shows that the distribution of sites is nearly Gaussian, implying 
that interactions between the polymer and dopant are small. 

Since broadening of the nonlinear susceptibility is shown 
to have an important affect the dispersion, the determination 
of excited state properties of molecules from spectroscopy 
requires that such a theory be used. So, IB theory using 
stretched Gaussian statistics may become an important tool 
for interpreting nonlinear-optical spectroscopy measurements. 
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excited state and the superscript * on W indicates the complex 
conjugate of ft. This allows us to use a simple frequency 
argument of ±u> that significantly improves the readability 
of the equations in the extensive appendices that follow. An 
example of the compact notation is given below: 



w 



(1) 



l2g 



(17) 



Since all arguments are all of the form 
-(fi» a T^.) an( j a jj n-ausjtiQjjs are from/to the ground state (g), 
this form describes all inhomogeneous broadening terms in 
this paper. Table [XT] summarizes the compact notation. 

TABLE II 

Compact form of wi x '(z) and W(z) up to third 

ORDER(x= {1,2,3}). 
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Appendix A 
Compact notation 

The energy denominators D for the higher-order susceptibil- 
ities are complex combinations of wi x \z) or W{z). We have 
developed a more compact notation than previously used lfTTl . 
For (3 < 1, we have added a subscript to (3 to indicate the 
excited state involved in the process and a * on the power of 
W to indicate a complex conjugate of the complex argument 
il. Similarly for (3=1, the subscript on W indicates the 
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Appendix B 

Energy Denominators for Selected Processes 
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i 



'/3m 



TO-") + W Z m} {< ) (-^) - < r m}' 

+ «J(-a;) + <f (2a;)} {w#>(0) - <>* ( W )}] } 



(36) 



8 



D/^ n (-2w;w,w,0) = 



^7T 



3/2 



r )lg^fmg'~ing 



{w^(-2u,)W^\-w) [w^(-2co) + w£>'(0) 
+W% r (w) w£f (2w) [w£ } (0) + W£ r (2w)" 

+W^*M [w£>(- W ) + w£>*(2u,)] [^ ) *(o;) + ^ 1) *(0)]} 



(37) 



2) (3= 1: 



2r 



7is \w(w + 2irj fl )7j 2 fl 
+■ 



2r 2 



w 2 (w + 2«r ;s ) 2 



[W,(-w) + - Wj*(0) - W,(0)] 



~o [Wl(-2w) + W,*(2w) - Wi(-w) - VK z *(a 



or 



(38) 



£>/,f(-2a;;a;,a;,0) 



— ^ [Wn(-w) {Wi(-2w) - W t (0)} + W*{u) {W?{0) - W 1 *(2uj)}} 

llgjng I 2UJ 

+ 1 [{W n (-u) + W n (0)} {W t (-2u>) - Wti-uj)} + {W*(u) + W:{0)} - W;(2u;)}} 

+ - ^ , [W(0) + W7(0)}{W„(- W ) - 



2(w + ir„ g ) 
1 

(w + 2iT ng ) 



[{W n (0) - W*(0) + W„(-o;) - W„>)} + W?(u;)}] 



(39) 



D/^(-2w;a;,w,0) - 



\W m {-w) {W,(-2w) - W0)} + W* M W(0) - W7(2w)}] 



+ - [{W m (-2u) + W m (-w)}{Wi(-2a;) - + (2w) + (w)} {W,* (w) - W7(2w)}] 



+ 



1 



2(w + ir /3 ) 

+ 7 . U ; [{W m (-2u) + W* m (a;)} {Wi(-w) - Wj*(0)} 



(w + 2ir (9 ) 



+ {W m (-w) + (2a;)} {Wj(0) - 



(40) 



£>fmn(-2w;W,W,0) 



r 3/2 



{W m (-2 W )W„(- W ) [W,(-2w) + W,*(0)] + (2w) [W„(0) + W*(2w)] 



+W ro (-u;) [W,(-2w) + [W„(-w) + W„(0)] 

+w* m [w n (- w ) + w n *(2 W )] + wno)}} 



(41) 



D. Kerr Effect 

1 ) Intensity dependent Refractive index: 
a) /3< 1: 



D{F(-u;u,-u,u) = 



Hg 



1 



7 2-<(--) 



+ 



?r 2 



W + 2zr; 9 

2iwr ; „ 



U) 



2iTi g ui 



4r 2 



is" 



wi 1} M + w^ (-w)-W^ (w) 



.(i) 



7*s 



(i)* 



7*3 



Pi 



(42) 



9 



n IB ( 

U ln V 



■u>: u>, u>, -co 



llglng 
1 



wj£(-w) + W<£{u>) 



+ 



2iT 



ng 



W${-w) + W£{u>)-W£ {-<*>)- W£ (co)\ [W^(-uo) + W£ H 
W^(-co) < r (co)] \w£\u,) + W$>\-w) 



-(1) 



llg y ' 



W 



n) '(u>) + w£>\-w) 



llg ^ ' 



(i)* 



,(1)* 



1 



2(co + iT ng ) 
1 



2co 



(43) 



— 7T 

llglmg 
1 



7i 9 



+ 



2(w + *r (9 ) 



l 

"2^ 



{w^(-2co) + <f(0)} {w£\-w) - W^(-co)} 
+ «f (2c) + « } M - < r M}] 

O) + <!>)] K )( - w) - < r H 

£ ) (°){<^)-< 1) m} + 



^ (2 W ) + W«>) 



7ig 



'(0){w^(-co)-W^ 



(44) 



Dii n (-uj;uj,uj,-Lo) 



-in 



3/2 



{O) + [< 1) (-c) + Ty«*M 



+W« (-2a;) Wj» (-co) [w£> (-co) + WV' (-co) 
+ W^'(2co)W^(co) [w£(co) + W^'(co)] } 



(45) 



b) = 1: 



Du B (-u;uj,-co 7 uj) = 



iV^ (2(n lg -w) 2 -r l 



Hg \ 



is 



+ 



+ 



21% - uo 2 + 2iTi g L0 



4i> 2 



\Wi(co) + w;(-uj) - w;(co) - Wi(-co)} 



oo + 2iYi g 
2icoTi„ 



(46) 



[W n (-w) + W n (w)] 



2 ( n y"W l( -") + — 



V^llg 



llglng 

+ [w:(uo) + w:(-uo) 

+ 2^~ i W n(~^ + W n(") - W:(-co) - W*(uj)] {W t (-co) + W;(co) 



V^llg 



- ng 

1 



2(co + iT ng ) 



[W n (-co) - W*(co)] [Wi(co) + W?(-u>)] 



1 



+ — [W n (-co) {W^-co) - W t (co)} + W*(co) {Wf(-co) - W t *(co)}} 



(47) 



10 



[W m {-2cu) + W m (0)} 



2(n, ff -w) 



Wi(-w) 



2i 



llglmg 

+ [W*(2 W ) + W*(0)] 

+ 2^ [{^-(- 2w ) + - ^'(-w)} + {^* (2w) + VK m (0)} - 



2(w + iT lg ) 



[W m (0) + W* m {Q)] [Wj(-w) - 



+ ^ [W m (0) {Wi(-w) - + W* (0) {WK-u) W*(u)}] 



(48) 



A J m n (-w;w,w,-w) 



{W ro (0) [W„(-w) + W n (w)] [Wi(-w) + 

ilgimging 

+W* (0) [W n (-u) + W*(u)] + W;(-uj)] 

+ W m (-2u)W n (-w) [Wti-uj) + + W^WM [W n (w) + W*(w)]} (49) 



2) Pump-Probe: 
a) < 1: 

Di l B (-ui;w 1 ,U2,-U2) = 



2i^/n ( 2uj x 

llg \ (Wl + W 2 )(Wl - W 2 ) 



_^)* (wi) __L VF (2) ( _ LJl) 



7i 3 



7/9 



+ 



+ 



(uji + u) 2 ) 2 {uji - lo 2 ) 2 {ui + ll>2 + 2iTi g )(uJi - lu 2 + 2iT ig ) 
(w§ - + MT lg w% - iYigUJibj 2 - 3r 2 ff u; 2 + tfgLJi) 



W$>{-<*) + W£> (wi) 



(i)* 



+ 



zr (g o) 2 (wi - w 2 ) 2 (wi + w 2 + 2«r ;9 ) 
+ - 3ir ;g w| - ir ;s u;iw 2 - 3r 2 g ^ 2 - r 2 g wi) 
zr ;s o) 2 (wi + ^ 2 ) 2 (^i - w 2 + 2ir (9 ) 



W^(-W2) + IV^ M 



W£\c02) + W£ r (-0J2) 



(50) 



An (-wi;wi,w 2 ,-w 2 ) = 



7/ff7ng 



7; 9 



+ 



+ 



+ 



+ 



1 



wi - w 2 + 2ir„ 9 

(1), 



W^(-W2) + W^(W2) 
W £\w2) + W£f{-U2) 
W^{-W2) + W£'{W2) 



^'(^O+^W-^ (-W 2 )-W^ (Wi) 



(1)* 



1 



<»( W2 )|< r (- W2 ) 



(1)*, 

01 



ui+u 2 + 2iT ng 



W^(-w 2 ) + W^(-wi) - (W2) - W^* (wi) 



1 



2^r 



ng 



<'U) + ^' (wi) 

{ - < ) (-^ 2 )] [< } (-^l) + Wj»(W2)] 

W^'^2) - < r ( Wl )] [<f ( Wl ) + <f (-c 2 )] } 

{ « 1J (- W l) - w#>m] [^(-Wl) + <^ 2 ) 



Wi — w 2 



+- 



1 



+ 



+ 



(i)* 



0i 



(_a*)-w£>*( W i)] [^'(wiJ-W^^wa)]}} (51) 



11 



CK; w i. w 2,-w2) 



llglmg 
-1 



TV 



- w 2 ) + Wp^(—wi + lo 2 ) 



1 



ui- bj 2 + 2iTi 



TV 



(i)" 



(wi + w 2 ) + (wi - w 2 ) 



{ [W<£' (0) + TV« (- Wl - lo 2 )] [lV« {-W!) - < r (-wa) 



Tvi 1} (0) + (wi+wa) 



(i)* 



< ) (c 2 )-TV«*(c 1 )]} 



+ 



1 

wi + w 2 + 2zV fl 

+ 



{[wV\o) + wj£(-w 1 +u 2 ) 



(-W2) 



+- 



Wi — w 2 



^ { [^(0) + ^(-wi - W2)] [w£\-ui) ~ w£\-uj 2 )] 



Wi 1} *(0) + Wi ir (a;i + 



+- 



+ [W^ (-wi + c^ 2 ) + TV £2* (wi + W2) 
1 



TV, 



Wl + w 2 



{ [<? (°) + <L J (""I + ^)] « 1} (-wi) - W$ (a*)] 
+ [wj»\o) + W<»' fa -<#>)] [w^'(-u 2 )-W^'fa)]}} (52) 



Amn(- W i; w l> w 2,-W2) = 



ft; /3 = 1: 

£»h B (-wi;o;i,w 2 , -w 2 ) 



-iTT 3 / 2 



{w^(-o;i - wa) [TV^(-a;i) + w£ r (-u; 2 )] [tV^-wi) + Tv£ } (-o; 2 ) 



+w£>(- Wl+W2 ) [TV«(- Wl ) + TV«*(a; 2 )] W'W+IV^^) 



+TV«(0) W^'i-uJ + Wf/ fa) TV^(-^ 2 ) + TV^V 2 ) 



+w 2 ) 

+Tv£ r (^i-^2) 



TV«*( Wl ) + TvW*(o; 2 ) 
< ir K) + < ir (-^ 2 ) 



TV«(o; 2 ) + Tv£ r ( Wl ; 



2^ J 



TV 



M + TV«*(-^ 2 ) 



W^(-u 2 ) + W^'fa) 
(-wiJ + W^Vi)]} 



(53) 



2wi 



2 < si '—) H .- I (-.)-^^v>7(^) 



+ 



+ 



+ 



7/9 I fa + w 2 )(wi - w 2 ) 
2 fa- 2 u> 2 w| + w| + 2ir 2 + 6r 2 w 2 - 8ir Jg wiw 2 ) 



Ti9 



(wi + w 2 ) 2 (wi - w 2 ) 2 (^i + lu 2 + 2iTi g )(u>i -uj 2 + 2iT lg ) 
fa - uiwl + 3ir ifl w| - iTi g uiu> 2 - 3r 2 9 w 2 + Tf g uji) 



[W,(-wi) + Wi*(wi)] 



ir /g w 2 (wi - w 2 ) 2 (wi + w 2 + 2«r ;s ) 
(w| + - 3zVgw! - ir /g wiw 2 - 3r 2 g ^ 2 - r 2 g ^i) 

?r ;g w 2 (wi + to 2 ) 2 fa - w 2 + 2iV s ) 



[TV ; (-co 2 ) + TV/M] 
[Wi(w 2 ) + Wr(-W2)]l (54) 



D I l ^{-u} 1 ;w 1 ,w 2 ,-u} 2 ) 



— IT 



llglng 



2(Qi g - mi) 



<4 



1 



wi - w 2 + 2iT ng 



w t *fa) + 



TV(-wi) + 
2i 



2i 



[TV„(-w 2 ) + TV„(w 2 ) 



V^Hg 



V^llg_ 

[W*fa) + W:{-u 2 ) 



lji + uj 2 + 2iT ng 



[W,(-W2) + [W„(-wi) + W n fa) - TV„*(-o> 2 ) - W*fa)\ 

[TVM + W,*(-W2)] [W n (-W2) + W n (-wi) - TV„*(^ 2 ) - W*(wi)] 



12 



+ 2lT~ + lWn(U2) + W n (-u 2 ) - W*{-w 2 ) - W*{oJ 2 )] 



+- 



- rig 
1 



UJi — U) 2 



+ - 



1 



Wi + w 2 



{[Wj(-Wl) - Wj(-W2)] [W„(-Wl) + W n (W2)] 

+ [W7(W2) - [W*(a;i) + W*(-W2)]} 

{[Wj(-Wl) - Wj(W2)] [W„(-Wl) + W„(-W2)] 



+ [Wl*(-W2) - W7(wi)] [W*(wi) - ^„> 2 )]}} 



(55) 



Di£ 1 (-uj 1 ;uj 1 ,W2,-uJ2) 



llglmg 



2(n tfl -^i) 



+ 



2(nr fl + ^i) 
1 



W,*(wi) + 



2z 



2i 



[W m (-wi - w 2 ) + W m (-wi + w 2 )] 



V^llg 



V^llg_ 



+ — {[W* m (0) + W ro (- Wl - W2)] - Wl*(-W2)] 

Wl — W 2 + ill Ig 

+ [W m (0) + W^{ui + wa)] [W,(W2) - W*(wi)]} 
+ ~ — 9 , r {[^* (0) + W m (-W! + lo 2 )} [Wi(-wi) - W7(a*)] 
+ [W ro (0) + W£(wi - wa)] W2) - W,*(wi)]} 

+ " {[Wm(0) + - W2)] [Wi(-Wi) - W,(-W2)] 

U^i — UJ2 

+ [W* m (0) + W* (wi + wa)] - 

+ 7T7T- {PC^i - W2) + - wa)] - W*(-^)] 

2«1 , 



+ - 



1 



UJl + u 2 



+ [W m (-wi + co 2 ) + W^cd + wa)] [Wj(w2) - W,*(W2)]} 
{[W ro (0) + + wa)] [Wj(-wi) - 



[WA(0) + W* (wi - wa)] [Wi*(-W2) - (o>i)]}} 



(56) 



A 7 mn(- W i; w l> w 2,-W2) = 



r 3/2 



{W m (-Wl - CC 2 ) pVj(-Wl) + W*(-W2)] [Wn(-Wl) + Wn(-W2)] 



+W m (-to 1 + u 2 ) [Wi(-wi) + Wf (wa)] [W„(-wi) + W„(W2)] 

+ ^ m (0) + [W„(-W2) + Wn(W2)] 

+ co 2 ) \Wt{^) + W;(cu 2 )} [W n {oj 2 ) + W*(wi)] 

+TC(W1 - W2) + W,*(-W2)] [PV„(-W2) + W,>i)] 

+W* m (0) [W?(u; 2 ) + Wi*(-wa)] [W n (-wi) + W^i)]} 



(57) 



